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Abstract 

Despite the fact that the problem of interference mechanism of magnetoresistance in semicon- 
ductors with hopping conductivity was widely discussed, most of existing studies were focused on 
the model of spinless electrons. This model can be justified only when all electron spins are frozen. 
However there is always an admixture of free spins in the semiconductor. 

This study presents the theory of interference contribution to magnetoresistance that explicitly 
includes effects of both frozen and free electron spins. We consider the cases of small and large 
number of scatterers in the hopping event. For the case of large number of scatterers the approach 
is used that takes into account the dispersion of the scatterer energies. We compare our results 
with existing experimental data. 
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I. INTRODUCTION 



At low temperatures the conductivity in semiconductors is supported by carriers hopping 
between locahzed states on the impurities. At low enough temperatures the characteristic act 
of hopping occurs not between the neighboring impurities but between the closest impurities 
in the thin energy strip around the Fermi level. This phenomenon is known as variable 
range hopping conductivity and is characterized by the well-known temperature dependence 
of resistance: 



where To is a constant and a is less than unity. When constant density of states g[E) is 
discussed, the exponent a is equal to + 1), where d is the system dimension. However 
when the conduction strip becomes narrower and the Coulomb gap becomes important, 
then a = 1/2 for any system dimension. The variable range hopping conductivity is, in 
particular, discussed in details in 

The problem of magnetoresistance in the semiconductors has been actual since 1950-th 
Q . Then it was understood that at strong magnetic fields the magnetoresistance is positive 
and can be exponentially strong. This magnetoresistance is due to electron wavefunction 
shrinkage in the magnetic field. 

However it has appeared that at weaker fields the magnetoresistance becomes negative. 
This phenomenon has been first described by Nguen, Shklovskii and Spivak [3] (for detailed 
review see It is related to the interference. During the hop the tunneling carrier suffers 
under-barrier scattering on the impurities that are outside of the conduction strip. So the 
resulting hopping amplitude becomes the sum of different tunneling paths, which interfere 
with each other. 

It is important that the hopping conductivity is controlled by percolation and, thus, by 
the largest resistors (formed by the pairs of hopping sites) that still should be included into 
the percolation network. This fact strongly emphasizes a role of destructive interference 
which can practically cut the percolation path and leads to an increase of the resistance. 
The magnetic field suppresses the interference and leads to negative magnetoresistance. 

The negative magnetoresistance appears to be linear (on the magnetic field) in the rel- 
atively wide range of weak fields, so it dominates over quadratic positive wave- shrinkage 
magnetoresistance at weak magnetic field. 




(1) 
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The theory of interference magnetoresistance [4] is based on the model of spinless elec- 
trons. This model can be justified when all electron spins are frozen, for example, by the 
exchange interaction. In it is stated that there is no interference effect on the magnetore- 
sistance when all electron spins are free. However this statement lacks the solid theoretical 
proof. Also, no discussion of the systems in which some spins are frozen and some are free 
is given. Therefore the question remains "how many free spins one should have to suppress 
the interference magnetoresistance?". 

In it is stated that the admixture of free spins leads to additional positive magnetore- 
sistance due to spin alignment in the magnetic field. However no detailed discussion of this 
magnetoresistance mechanism is given. For example its detailed temperature dependence 
remained unknown. 

In the studies 3|, |J] authors consider the system with a large number of scatterers involved 
in a hopping event. In our study we denote this situation as the case of long hops. For this 
model only qualitative results are obtained in [3|, IJ]. Later it was shown that in 

many realistic systems the characteristic number of scatterers is small. Most hopping events 
occur without scattering, and to consider interference magnetoresistance one can take into 
account only hopping events including a single scatterer. We denote this situation as the 
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case of short hops. For this problem of short hops (and spinless electrons) the papers 
give quantitative estimates of the magnetoresistance. 

There were several attempts js iS] to develop a quantitative theory for long hops. How- 
ever, to our opinion, no one of this works pays a sufficient attention to the variance of scatter 
energies. 

Also we state that the conventional theory of interference contribution to magnetoresis- 
tance is not sufficient to describe existing experimental data. While semiconductors with 
hopping conductivity usually exhibit a combination of negative linear and positive quadratic 
magnetoresistance (this fact is in accordance with conventional magnetoresistance theory), 
the temperature dependence of negative and positive magnetoresistance differs sometimes 



from predicted by existing theories 13l-ll9| . In 13 
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spin effects. In 
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191 this difference is attributed to 
19 1 the spin alignment magnetoresistance (suggested in j^) is invoked to 



describe experimental observations. However no rigorous theory of this magnetoresistance 
(and any consistent discussion of free spins effect on the magnetoresistance) is presented. 
The goal of the present study is to develop theoretical approach to the interference con- 
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tribution to magnetoresistance that explicitly takes into account an admixture of free spins. 
We consider two general problems. The first is to what extent the interference magnetore- 
sistance is suppressed when spins become free. The second is related to derivation of the 
expression for spin alignment magnetoresistance that would give the detailed dependence of 
this magnetoresistance on the system parameters including magnetic field, temperature and 
the weight of the free spins. 

In our work we consider both long and short hops. For long hops (large scatterer numbers) 
we use an approach that takes into account the scatterer energy variance. Our approach 
becomes rigorous for exponentially broad distribution of scatterers energies. However we 
show with numerical computations that it is applicable to realistic energy variances. 

The plan of this study is as follows. In section [ITl we discuss the basics of interference 
magnetoresistance mechanism. In section IIIII we present theoretical analysis for the short 
hop case. In section [IV] we discuss the statistics of tunnehng path amplitudes in the long 
hop limit. In section |V] we use this statistics to get expression for magnetoresistance in this 
case. In section |Vl] we compare our results with existing experimental data and in section 
IVIII we present the final discussion of our results. 

II. INTERFERENCE MECHANISM OF MAGNETORESISTANCE 

It is convenient to describe hopping conductivity in terms of random resistance network 
first proposed by Miller and Abrahams [2^. In terms of this model any pair of impurities 
(denoting as 1 and 2) corresponds to a resistor with resistivity 



Here is the hopping rate, ni and n2 are the occupation numbers of the first and the 
second impurity states, correspondingly. /12 oc exp(— ri2/a) is the energy overlap integral 
between states 1 and 2, ri2 is the distance between impurities and a is the localization 
length. Nph is the probability to find a phonon for a hop (one can take Nph = 1 for hops 
with emission of a phonon). Finally, angle brackets mean time averaging. 

The terms rii, (1—722) and Nph specify the dependence of resistor R12 on impurity energies. 
We will not discuss this part of (|2]). What we are interested in is the overlap integral J12. 




(2) 
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When there are no impurities beside 1 and 2, the overlap integral can be easily evaluated [l| 



'12 



a J 



exp(-ri2/a). (3) 



Here Jq is the constant of the order of the Bohr energy, (5 is the pre-exponential factor that 
is determined by the impurity type. For shallow impurities in the bulk system (or in the 
5-doped layer) /3 = 1. For non-Coulomb impurities (like the deep impurities or doubly- 
occupied 01 centers) /3 = — 1. In two-dimensional systems /3 is equal to 1/2 and 
— 1/2 correspondingly for Coulomb and non-Coulomb impurities. 

When there are intermediate scattering impurities between the hopping ones, the overlap 
integral becomes a sum of the different tunneling paths and expression ([3]) stands only for 
one of them (the one without scattering). The total overlap integral between impurities 1 
and 2 can be expressed as 

/r. = /. + E^ + i:^ + ... (4) 

i * ij ^ J 

Here lij are defined by (j3]). Ei are the energy positions of scattering impurities (with respect 
to the Fermi level). We assumed that the corresponding energies are much larger then the 
energies of hopping impurities 1 and 2. Each summand in (jlj) corresponds to an amplitude 
of some tunneling path. Let us note that any tunneling path amplitude is oc exp(— //a) 
where / is the length of the corresponding hopping segment. So the paths with large I can 
be neglected. Accordingly, we neglect all paths with backscattering and also include in our 
consideration only scattering impurities in some thin strip (in real space) around the line, 
connecting impurities 1 and 2. The width of this strip is estimated in \^ as ^r\2a. However 
we will show that in the case of long hops it is sufficient to consider a thinner strip. So we 
have a finite number of paths 2^ where iV is the number of scatterers in the strip. 

Let us also write the expression for the simplest (but important) case of one scattering 
impurity s. In this case 

oc |Ji + J2|% Jl = /l2, ■J2 = ^^^- (5) 

Es 

Note that the expression (|5]) and the upper expressions do not take in account the electron 
spin. Consequently up to this moment we discussed the model of spinless electrons. As 
we have already mentioned, this model can be justified when electron spins on scattering 
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FIG. 1: Final results of direct and scattered tunneling path in different spin configuration. 
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impurities are frozen (for details see [^). But now let us assume that the intermediate 
impurity is occupied and its electron spin is free. 

The scattering on the occupied impurity is equivalent to the cotunneling of the initial 
electron to the scatterer and the electron from the scatterer to the destination impurity 
(figured]). One can see that when spins of initial electron and electron on the scatterer have 
the same directions, the final result of direct tunneling Ji and the tunneling with scattering 
J2 are the same. However, when the directions of these spins are antiparallel, the final 
results of the trajectories Ji and J2 are different. So in this case the tunneling paths do not 
interfere. 

As long as the spin on the scattering impurity is free, there is always a probability that the 
two spins have the same direction and a probability that they have antiparallel directions. 
Note that the expression ([2]) for Tij contains time averaging. So for the free spin on the 
intermediate impurity one obtains 

oc Pn 1^1 + -hf + (1 - Pn) (l^il' + l-hf) ■ (6) 

Here P^-^ is the probability for two spins to have the same direction. Without the magnetic 
field P-|-f = 1/2 because both spin directions are equally probable. The external magnetic 
field aligns the spins leading to an increase of P-^^. The expression for probability P^^ in a 
magnetic field is as follows 

Here fi is the magnetic moment of a localized electron, H is the external magnetic field and 
T is temperature expressed in energy units. 

In the expression ([6]) we have implicitly assumed that the spin states change with time 
faster than the hopping events occur or at least faster than Miller-Abrahams network is 
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stabilized. To explain this assumption we note that the spin states at least can not change 
significantly slower then the hopping rate because the hopping in the antiparallel spin com- 
bination can exchange the spins. Also, other mechanisms of spin diffusion or relaxation can 
exist, for example, related to the exchange of spin directions between neighbor impurities 
with free spin. We also note that our expression (jS]) agrees with expression (14) from jsl 
(that describes spin alignment magnetoresistance). 

There is another assumption that we will rely upon. Namely, we will make use of the 



fact that the frozen spins form the so called Bhatt-Lee phase [2l|]. In this state the frozen 
spins are in a singlet state and thus have no preferred direction. So the frozen spins scatter 
all tunneling electrons in the same way (independently of their spin direction). If by some 
reason the frozen spins form spin glass state, our considerations can not be applied. 

One can see that there is a dependence of on magnetic field even if the intermediate 
impurity has a free spin. First, in magnetic field a phase difference between tunneling paths 
exist. It can be expressed as J2 — )■ ^26*"^, ^ = HS/^o, where S is the area between the 
tunneling paths, $0 is the magnetic fiux quantum. Secondly, the probability P-^-f depends 
on the magnetic field (see Eq. ([71)). According to the logarithmic averaging procedure (j^j), 
one can write the following expression for the contribution of scattering on free spins states 
to the magnetoresistance: 

In ^^^^ - An '-^^ + + ~ ^^^^ ^''^^1' + '-^^''^ \ (8) 

^(0) \ \J^ + J2f/2 + i\J,\' + \J,\')/2 /^_- 

Here the angle brackets with index "free" mean the averaging over critical resistors with 

free spin on the scatterer (note that in short hop limit we have only one scatterer) . There is 

also a contribution from resistors with unoccupied intermediate impurity and from resistors 

involving a scattering impurity with frozen electron spin. These contributions have a similar 

form: 



R(0) \ \Ji + M 



frozen 



Note that angle brackets with index "frozen" mean not only the averaging over scatterers 
with an electron with frozen spin but also over the unoccupied scatterers. The contribution 
of these two scatterer types have the same form although the effects of these scatter types 
are different due to the sign of energy Eg. 

Now we generalize the discussed approach for the situation when a resistor includes N 
scattering impurities and M of them are occupied by electrons with free spins. When a resis- 



tor contains one free spin, there are four possible initial spin configurations (corresponding 
to up and down projection of initial spin and free spin on the scatterer). Similarly, M free 
spins in a resistor lead to 2^+^ possible initial spin configurations. Also, in a resistor with 
one free spin one has no more then two possible tunneling results (final states of the tun- 
neling) for any initial configuration — with and without exchange of spin projections. In 
our resistor one has no more than 2^ tunneling results for any spin configuration. These 
results can be described by the combination of free spins that changes the spin projection 
during the process of hopping. Note that some results are impossible in specific spin con- 
figurations (for example for one intermediate free spin there is only one tunneling result for 
spin configurations with parallel spin projections). 

Accordingly, for any spin configuration 2^ tunneling amplitudes can be separated into 2^ 
groups corresponding to different tunneling results (some groups can contain no tunneling 
amplitudes in some initial spin configurations). The tunneling paths that are inside one 
group interfere (because they lead to the same final state in the given initial configuration 
of free spins). Other tunneling paths do not interfere. One should average over initial 
configurations: 

2M+1 2 

(10) 



2M+1 

oc y: Pc{h) e 



c=l 



k{c,r) 



Here index c enumerates possible spin configurations, Pc{H) is the probability of a given 
spin configuration in the magnetic field to exist. Index r enumerates the tunneling results. 
Jfe are the amplitudes of the tunneling paths and (pk are the corresponding phases. The sum 
over k is taken only over those tunneling paths that lead to result r in the spin configuration 
c. The probability Pc{H) can be expressed as 

exp [{2N^p{c) - M - 1)^H/T] 



Pc(H) = 



(11) 



2M+icosh^+\nH/T) 

where Nup{c) is the number of free electron spins (including the spin of tunneling electron) 
that are aligned along the magnetic field in the configuration c. 
The corresponding expression for magnetoresistance is 

2 



RjH) 
R{0) 



In 



2M+1 

E Pc{H) E 

c=l r 



fe(c,r) 



2M+1 

E PM E 



c=l 



E 

k{c,r) 



(12) 
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Here the angle brackets mean averaging over critical resistors. We note that number M, 
possible spin configurations and tunneling results can be different for different critical resis- 
tors. 

When there are no resistors with two or more scatterers, the expression f[T2]) is reduced to 
(jH]) and (jH]). When all electron spins are frozen, it is reduced to the conventional expression 
for magnetoresistance from (that corresponds to oc |/i2p with J12 defined in (jl])). 

The expression f|T2|) is rather complex. However it allows numerical computation of mag- 
netoresistance for a relatively large number of scatterers. The results of such computations 
will be presented in the section |V] of this study along with analytical analysis of the magne- 
toresistance. 



III. SHORT HOPS 



In the case of short hops one can neglect resistors with more than one scatterer. Hopping 
events with no scatterer do not contribute to interference magnetoresistance so that one can 
consider only resistors with one scatterer. For those resistors it is possible to use expressions 
(ED and (P. 

The magnetoresistance corresponding to short hop limit with no free electron spins (ex- 
pression ([2])) is known from [g], Q]- It can be shown that the averaging in yields dif- 
ferent results for the two areas. In the first area, Ai, the denominator is not very small 
I Jl + J2I > fJi- In this area one can expand the logarithm. It leads to a quadratic magne- 
toresistance. In the second area, A2, where | Ji + J2I < y^Ji, the logarithm is large and can 
not be expanded. This area leads to the negative linear magnetoresistance that dominates 
in weak fields. For short hops the expression for this magnetoresistance is (see 6|, [7|) 

R{H) 2+d/2Tj /1Q\ 

m — -— oc — r;, ' H, (13 

where Vh is the mean hopping distance and d is the system dimensionality. The dependence 
on Th controls the temperature dependence of the magnetoresistance as rh oc T^^^^'^^^^ for 
the Mott-like hopping and Vh oc T~^/^ for Efros-Shklovskii hopping over Coulomb gap states. 

For the free electron spins the area A2 does not exist. Actually, the denominator in ([8]) 
is not less than 3|Jip/4. The denominator in (Q for resistors with unoccupied scatterers 
(that has constructive interference) is larger then Hence in weak magnetic field the 
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logarithm can always be expanded. This leads to the following expression for the weak field 
magnetoresistance expansion: 

R{H) ^ I {^fJ,.h - \-h-h^' \ lJ^hf_\ 

Here angle brackets with index "occupied" mean averaging over resistors with occupied 
scattering impurity (in this expression we consider all such impurities to have a free electron 
spin). Note that for this impurities J1J2 is negative. Angle brackets with index "unoccupied" 
corresponds to averaging over resistors with free intermediate impurity. For those resistors 
J1J2 is positive. 

There are three contributions to magnetoresistance due to interference for free electron 
spins. First, there is one negative contribution. It comes from magnetic field suppression 
of destructive interference for the resistors with occupied scatterers. However, in contrast 
to the situation when electron spins are frozen, it is quadratic (not linear) in terms of the 
magnetic field and does not automatically dominate over other terms in a weak field limit. 
Also, it does not automatically dominate over the wave shrinkage magnetoresistance even 
for weak fields. 

The second term appears from the suppression of interference in resistors with unoccu- 
pied scatterers and constructive interference. It will dominate the first term at least for 
semiconductors with large compensation K {1 — K 1), a.s the number of unoccupied 
scatterers in this semiconductors is larger than the number of occupied ones. 

The third term in ( |T4|) is also positive and quadratic in terms of magnetic field 

Here n is the impurity concentration, d is the system dimensionality. On the right-hand side 
of the equation we used the following approximation. We considered Ji and J2 to be of the 
same order (it is usually correct if the scatterer lies in the thin area between the hopping 
impurities |^). the probability to find a scatterer in this area. 

The term (03]) is related to the dependence of the probability P-j^^ on the magnetic field. 
It is the spin-ordering magnetoresistance first discussed in jsl. One can see that this term 
has very strong temperature dependence. For example, for a 2D system with hopping over 
the Coulomb gap states, this term depends on the temperature as T~^^/^. To compare, the 
temperature dependence of the wave-shrinkage magnetoresistance in this case is oc T~^/^. 
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Therefore at low temperature the spin-ahgnment magnetoresistance should become stronger 
than the wave-shrinkage one. 

However, at low temperatures the system should correspond to the case of long hops 
(rather then of the short hops) and at some temperature electron spins should be controlled 
by the exchange interaction. However (as we show in section |V]) the spin alignment magne- 
toresistance remains in the long hop limit and has a strong temperature dependence. 

Finally, one can consider the situation when some spins are frozen and some spins are 
free, that is < P/ree < 1- The short hops theory of magnetoresistance can be easily 
generalized for this case. Naturally, each resistor contains either a free spin or a frozen one 
(if it involves an occupied scatterer). The general expression for magnetoresistance contains 
averaging over critical resistors. Accordingly the following equation can be used for shot 
hops 

In ^(^Q^ ~ PfreeMRfree + (1 " PfreejM Rf rozen, 

where MRfree is given by the expression f|T4l) and MRf rozen is given by the expression (fT3ll . 
Naturally, when there is a significant number of frozen spins in the system (1 — Pfree is not 
very small) the weak field magnetoresistance will be controlled by MRfrozen and will be 
negative and linear in terms of magnetic field. However, if one considers the quadratic term 
of the magnetoresistance expansion, then the term MRfree can also become important. 



IV. LONG HOPS. STATISTICS OF TUNNELING PATHS 

We address now the case of long hops, i.e. the situation when there is a large number 
of scatterers in a hopping process. As we have mentioned before, in this case there is 
exponentially large number of tunneling paths 2^ (even without backscattering) , where N 
is the number of intermediate impurities in the resistor. When there are no free spins, all 
paths interfere. The situation with free spins is more complex as it was discussed in section 

m 

In our opinion, the role of free spins in this case is even more important than it is for 
short hops. Indeed, for sufficiently long hops free spins are always present in the area 
between hopping sites. It is known that with decreasing temperature the spins of electrons 



that are localized on the impurities form the so-called Bhatt-Lee phase 21|. One of the 



properties of this phase is a relatively large number of free spins. The simplest estimate of 
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the number of free spins is based on the nearest-available- neighbor approach 22|. It leads 
to the logarithmic decrease of the free spin concentration with decreasing temperature. The 
hopping length grows with temperature according to a power law. So at low temperatures 
the number of scatterers with free spins in a hopping process is large. 

In {4] the following arguments for the absence of interference magnetoresistance in the 
case of free electron spins has been proposed. It has been stated that in this case there are 
many tunneling paths that do not interfere with each other, making the strong destructive 
interference (that is responsible for linear negative magnetoresistance) very unlikely. 

However, at sufficiently low temperatures each critical resistor has free spins and, there- 
fore, there is always a large number of tunneling paths that do not interfere (there are 2^ 
possible tunneling results, where M is the number of scatterers with free spin in the hopping 
act). Consequently the quantitative estimate of the effect of the free spins on interference 
magnetoresistance required. The present study provides such an estimate. 

The role of interference effects in the magnetoresistance in the limit of long hops has 
been discussed (although without free spin consideration) in a number of publications 



121]. However all these reports imply two important model assumptions which we do not 



believe to be obvious. 

The first assumption is placing of scattering impurities in the nodes of some lattice. 
Usually one places scatterers in the nodes of a finite square lattice. The hopping impurities 
are on the diagonally opposite corners of the lattice. Then, one considers the shortest 
hopping paths that are along the lattice bonds (there are many such paths with equal length). 
In the present study we do not discuss the applicability of this assumption. However we do 
not use it. 

The second assumption implies a neglect the variance of scatterer energies. In 

m 
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12| 



the authors consider scatterers with two possible energies W and —W, so that the variance 



of the absolute value of the scatterer energy was neglected. In |8|, 111] the authors consider 
the fiat distribution of scatterer energies at least as one of the possible models. However, 
then they assumed that the distribution of tunneling path amplitudes is the Gaussian one 
with the relative variance of the order of unity. 

We argue that the second assumption is not correct for realistic semiconductors. Actually, 
the scatterer energies in semiconductors have fiat distribution and this fact leads to the log- 
normal distribution of tunneling path amplitudes which is significantly different from a 
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normal one. This distribution plays a crucial role in our theory. 

Let us consider some tunneling path. In accordance to the expression (jl]), its amplitude 
can be given as 

= n (16) 

k{i) 

Here Ji is the amplitude of i-th tunneling path. The index k{i) enumerates the scattering 
impurities that participate in the tunneling path i. A^j is a number of these impurities. 
Usually Ni ~ N/2, where is the total number of scatterers in the resistor, ki is the first 
of the impurities of this path. Index (/c = 0) stands for the initial hopping impurity. The 
last index k = Ni + 1 corresponds to the destination impurity of the hop. 

Let us take a logarithm of the absolute value of J, and write explicitly the term that is 
related to the hopping distance. 

In -= -r;,/a + ^/i^, fil = + ln— - + /31n (17) 

-'0 \^k\ fk-l^k+lO' 

Here Ar^ is the additional distance an electron should tunnel due to scattering on the 
impurity k, Ar^ = + rfc fc+i — rfc_i ^+1. is the distance between scattering impurities 
i and j, it may occur that some of this indices {i and j) are equal to zero or to A^j + 1 which 
means the starting impurity of the hop and the destination impurity, correspondingly. (3 is 
the pre-exponential factor in the overlap integral (see ([3])). 

We consider the scattering impurities to be situated in a thin strip to prevent the addi- 
tional distance to be too large, Ar^ < a. Also, the typical values of scatterer energies are of 
the order of the Bohr energy, so Io/\Ei^\ ~ 1 (the sign of is random). Finally, the third 
term is of the order of In r^c/ a where r^^ is the characteristic distance between scatterers in 
the resistor. In real situations this logarithm is not very large. Accordingly, we can consider 
jj^f^ as a random value with expectation and variance of the order of unity. 

Sometimes it is more useful to sum not over scatterers participating in the path i, but 
over all scatterers in the resistor. Therefore we introduce the modified values /I* that are 
equal to /i* if scatterer j participates in the given path and equals to otherwise. 

(fi), j G path i ^ ^ 

(18) 
otherwise 

Again, we consider as random quantities with math expectation and variance of 
the order of unity. 
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If one takes /I* as independent quantities (this assumption will be discussed later), one 
gets the normal distribution for ln|Jj//o| (in accordance with the central limit theorem). 
The math expectation of In | Jj//o| is —r^ja + A^-E^ and its variance is ND^. 

Hence the quantity | Jj| has a log-normal distribution law 



(ln(|J,|//o) +Va-iVE, 



\2 



The math expectation E{\Ji\) and variance D{\Ji\) are expressed as follows: 



Jil > 0. (19) 



E{\J,\) = Jo exp [-^ + NE^ + ; (20) 

D{\J,\) = I^exp (^-^ + 2NE^ + 2iVD^) . (21) 

This distribution law is significantly different from the one that can be obtained from 
normal distribution of J^. Actually, it is the exponentially broad distribution. The important 
feature of the law ( fT9l) is that its math expectation (|20l) is exponentially larger than the value 
corresponding to the maximum of distribution density /q exp (— r/^/a + A^-E^). It means 
that a sum of many tunneling path amplitudes (for example /12) should be dominated by 
exponentially small number of tunneling paths. In the next section we use this fact to derive 
a theory of interference magnetoresistance in the limit of long hops. 

However, before we will discuss the magnetoresistance, let us consider the following mo- 
ment. We derived the log-normal distribution ( |T9|) from the assumption that the quantities 
/I* where independent. Strictly speaking, this assumption is incorrect. The same impu- 
rities that form the interference pattern of the resistor participate in all tunneling paths 
(although in different combinations). Consequently we believe that it is reasonable to check 
the log-normal distribution by numerical calculations. 

Let us take the resistor with 10 scattering impurities. This resistor has 2^° = 1024 
tunneling paths, which is enough to reach the reliable statistics. The corresponding results 
for tunneling path amplitude distribution are shown on fig. O One can see that they are in 
a good agreement with log-normal law. 

Another important property of the obtained tunneling path amplitude distribution is the 
fact that the sum of all tunneling paths is dominated by a small number of the summands. 
Let us in addition verify this result. In order to do this, we have to formalize the concept 
" dominate" . We do this in the following way. 
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FIG. 2: Distribution functions of Jj (a) and In Ji/Iq (b) calculated for a tunneling path amplitudes 
in a single resistor compared to the log-normal distribution. 

Consider the sum of all absolute values of tunneling path amplitudes J2i \Ji\- This sum 
has 2^ summands. Let us now take some number Assign of its largest summands, so that 



E l-^fel > 0.6^1 J.I 



(22) 



The left hand part of the inequality is the sum of Nsign largest | Jj | . The Nsign is considered 
to be the smallest number to satisfy inequahty ( l22l) . We will specify Ngign tunneling paths 
with largest absolute values of amplitude as the significant paths. 
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FIG. 3: (a) — dependence of Ngign on the scatterer number N. (b) — relative part of the 
significant paths Nsign/'^^ versus scatterers number. 



If the most of the tunneling paths are important within the sum X^i \Ji\i should get 
Nsign ~ 0.5 ■ 2^. However if our assumption is correct and X^i \Ji\ is controlled by a small 
number of its summands, we should have Assign ^ 2^. The corresponding numerical results 
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are shown on figure ([3]). One can see that, although the number Ngign exponentially grows 
with iV, the relative part of significant paths exponentially decreases with iV, demonstrating 
the validity of our assumption that the sum of contributions of different paths is controlled by 
a small number of summands. In the next section we will exploit it do derive the expression 
for magnetoresistance in the long hop limit. 



V. MAGNETORESISTANCE IN THE LIMIT OF LONG HOPS 

In the spinless electron model the interference magnetoresistance is controlled by the 
overlap integral J12, which is the sum of many tunneling path amplitudes Jj. We have 
shown that this sum is controlled by a relatively small number of the largest Jj. Now we 
consider the question which of Jj occurs to be large. Jj is proportional to the product 
Hj exp(/i*) for impurities j participating in path i. The values /i* control whether the path 
i gives larger or smaller amplitude with an inclusion of the impurity j. 

Strictly speaking, /i* is controlled not only by the properties of impurity j but also by 
other impurities participating in path i. However, we can estimate /i* for any scatterer j 
and a characteristic tunneling path. 

7^=-:r^ + ln^ + /31n^. (23) 

Here Uj is the distance between impurity j and the line connecting the starting and the 
final impurities of the hop, r^c is the characteristic distance between the scatterers in the 
tunneling path, r^c will be estimated later in our paper. 

There are three opportunities for /ly. First it may occur that JTj is negative and has a 
large modulus, so exp/ZJ ^ 1. The corresponding impurities usually decrease the tunneling 
amplitude. The inclusion of a large number of such impurities makes the tunneling path 
amplitude exponentially small, so that this path will not be significant. Let us call such 
impurities with small Jlj as the irrelevant scatterers. Actually all impurities that are far 
from the line connecting the hopping impurities are irrelevant. 

In contrast, it can occur that the absolute value \Ej \ is much smaller then Jq and expjij ^ 
1. Inclusion of such impurities makes the path amplitude larger. If some tunneling path 
does not include many such impurities, its amplitude becomes exponentially smaller than 
amplitudes of other paths (which include the corresponding impurities) and this path can 
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not be significant. So the impurities mentioned above should exist in most significant paths. 
We call these impurities as the backbone impurities. 

Finally, there are impurities with exp/J^ ~ 1. They can be included and can be excluded 
from a significant path. We call these impurities as the interference impurities. 




FIG. 4: Dependance of mean numbers of the different impurity types (A'^;,;,, Nint and Nirr) on the 
total number of scatterers A^. Each dot is averaged over 100 realizations. 
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FIG. 5: Dependance of the relative number of backbone impurities Pf,;, on the system parameters, 
(a) — dependance on the width of impurity band AEhand with na^ = 0.2 and pre-exponent factor 
/3 = 1, (b) — dependance on the impurity concentration with AEi,and = lo and /? = 1, (c) — 
dependance on the impurity concentration with AEfjand = ^0 and /3 = — 1. 

Strictly speaking, the separation of impurities to irrelevant, backbone and interference 
ones is valid only for the exponentially broad distribution of scattering impurity energies 
(in this case Inlo/Ej is the leading term in ( l23l) and can be large). However, even if 
the distribution is not exponentially broad, there are some impurities with small value of 
\Ej\ that are good candidates to be the the backbone impurities. Let us try to find such 
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impurities in numerical computations. To perform the computations, we consider impurity 
to be a backbone one when 90% of the significant paths include this impurity. Similarly, 
the irrelevant impurities are the ones that are not included in 90% of significant tunneling 
paths. 

Figure m shows the averaged results for resistors with different length and with a constant 
width of the area occupied by the scatterers. One can see that the number of backbone 
impurities Nhb depends linearly on the scatterer number N. So we can estimate the relative 
number of backbone impurities in the long resistor limit. The same procedure can be followed 
for numbers of interference and irrelevant impurities N^nt and Nirr- Note that the relative 
parts of different types of impurities can depend on the width of the area exploited for 
calculations. For very wide area most impurities will be irrelevant. So we redefine the 
relative number of backbone impurities to make it model-independent 

P,, = lim (24) 

Pbb is the relative part of backbone impurities with irrelevant impurities excluded. It does 
not depend on the the width of the model area if it is large enough. Fig. |5] shows our 
computation of the dependence of Pbb on the impurity band width AEband, concentration n 
and pre-exponent factor /3. One can see that Pbb increases with the decrease of AEband- For 
large band width Pbb should tend to some finite value as the impurities with large energy are 
always irrelevant scatterers. Pbb weakly depend on concentration and pre-exponent factor 
at least in observed range of parameters. 

Now let us make use of the existence of backbone impurities to calculate the magnetore- 
sistance in the long hops regime. In the present study we restrict ourselves with the simplest 
model of independent links. This model assumes that any two interference impurities are 
always separated with a backbone one. Therefore this model should be valid for Pbb > 0.5. 
However the results of our computations agrees with this model even for smaller Pbb- 

Let us illustrate our considerations with the following model resistor (see fig. [6]). We 
consider, first, the spinless electron model. This resistor contains two backbone impurities, 
two interference impurities and two irrelevant ones. So the total number of scatterers is 

= 6, the number of tunneling paths is 2^ = 64. However, in our model we assume 
that one can restrict himself with consideration of significant paths — ones that contain all 
backbone scatterers and bypass all irrelevant impurities. Therefore instead of 64 tunneling 
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FIG. 6: The example of the hopping resistor with six scatterers. The scattering impurities are 
numerated while the starting and final impurities in this scatterer are marked with indexes 5 and 
F correspondingly. 

paths one can consider only four ones: (S"— J-S— S-S— i-F), (S— 5'2— i-S— 7-5— t-F), 
(5 ^ 3 ^ 5 ^ 6 ^ F) and (5 ^ 2 ^ 3 ^ 5 ^ 6 ^ F). 

The amplitudes of these paths are expressed by Eq. (ITTI) . One can see that all amplitudes 
of these paths contain the term /35/F3F5. It corresponds to backbone impurities. Then, all 
paths contain connections between the backbone parts. For example, impurities S and 3 can 
be connected with either 1 33 or 153/23 /F2. One can show that within our approximations 
these connections are independent, i.e. the net tunneling amplitude J = can be 

expanded as 



The backbone impurities effectively separate the resistor into independent parts — the 
interference links. The amplitude of each link can be described in the same way as a 
tunneling amplitude for the resistor with one scatterer. The net tunneling amplitude of the 
resistor in the independent links model is the product of terms corresponding to all links 
and additional terms corresponding to backbone impurities. 



Here the first term is the product of 1/ Ej for all backbone impurities j and the product of all 
the overlap integrals Ikj between neighbouring backbone impurities k and j (there should be 
no interference link between these impurities). The second product is over interference links. 
We assume that each link i contains starting and final impurities il and i3, correspondingly, 




(25) 





(26) 



backbone 3 interference 
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(that are the backbone impurities) and intermediate interference impurity i2. For example 
in our model resistor for the link connecting impurities 5* and 3 the index il corresponds 
to the starting impurity S. The impurities i2 and ^3 are actually the scatterers 2 and 3 
correspondingly. The expression fl2^ is the generalization of fl2^ for the arbitrary resistor. 

From Eq. f l26p one can easily obtain the expression for magnetoresistance corresponding 
to long hops and spinless electrons in our model of independent links. 

_ \MH) + MH)f\ ^ L \MH) + _ 



int 



where 



MH) = hiAH), j^, = IlldE^l^, (28) 



Here index i numerates interference links; Iii^i2, Iii,iz and Ii2,iz depend on the magnetic field 
through the phase (p. N^t on the right hand side is the mean number of interference links 
between the hopping impurities. It is proportional to the characteristic scatterer number 
A^. Angle brackets with index int mean that the averaging is over interference links (that 
correspond to interference scatterers in a resistor). 

is controlled by the area which does not give a significant addition to the tunneling 
exponent. Due to the presence of backbone impurities one should take the area with the 
constant width (independent on hoping distance) p ~ ^/r^ ^ y/r^a. The width p and the 
distance between neighbor impurities can be expressed as 

P~ - , r,,~ . (29) 



n/ \ n ^ 

Accordingly, the mean number of scatterers A^ = nrhp'^~^ is proportional to the hopping 
length Th. 

Finally, let us note that the quantity under the logarithm in (I27|) can significantly differ 
from unity only when amplitudes Jn and Jj2 are comparable (otherwise one of the amplitudes 
dominates and the discussed quantity is ~ 1). So we can substitute the averaging over 
interference links by averaging over all scatterers and, correspondingly, substitute Nint by 
N. 

m \ ihm + -Mo)f I ^ ' 

Here averaging is over all scattering impurities i. In the explicit expressions for Jn and Ji2 
( 128|1 one should substitute impurity i2 by impurity i. The impurities il and i3 should be 
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substituted with the backbone impurities neighboring to i. With this expression the problem 
of magnetoresistance in the case of long hops is deduced to the magnetoresistance in the 
case of short hops for resistors with ~ ^r^c- 

The expression fl5Ul) comes from the approximation of the sum of all tunneling paths J2 Ji 
and thus does not include free spins. To include free spins, one should derive the expression 
for magnetoresistance from (fT2!) with the approximation of A*";,;, ^ A^mt- The corresponding 
expression for is derived in the Appendix |X] (eq. IA3I) . With above mentioned arguments 
one can get from (lASp the following expression for magnetoresistance for long hops that 
includes the contribution of free spins: 



In 



:i-^..e)iV^#M±Ml^\- (31) 



R{H) 

m " V" \Mo) + Mo)\' 

_p ^ L Pn \MH) + MH)\' + (1 - ^tt) {\MH)\' + \MH)t 
\ \Mo) + .uo)f/^ + {\Mo)f + \Mo)f)/2 I 

One can see that terms in angle brackets in ( l3TI) have the same form as the magnetoresis- 
tance in the case of short hops (compare with expr. ([8]) and (IHl)). This magnetoresistance is 
discussed in section [XTTl However the "length" of the "resistors" corresponding to this terms 
in fl3Tl) is controlled by the distance r^c between neighbor scatterers in a resistor rather than 
by the actual hopping distance Vh- Note that Vgc is still larger than the distance between 
neighboring impurities n~^/'^. 

The expression (!3T1) means that all properties discussed for the short hops in section IIIII 
exist also in the case of long hops (at least as long as the independent link approximation 
remains applicable). There is a negative linear magnetoresistance due to the suppression of 
the interference by the magnetic field 

In ^ oc -iV(l - Pfree)rT"H oc -rr^H. (32) 

One can see that the dependence of the magnetoresistance on the hopping distance r/j (and, 
correspondingly, on temperature) becomes weaker than for the short hop case (compare 
Th to r^}^'^^'^)- At large temperatures there is additional temperature dependence due to a 
decrease of Pfree with decreasing temperature. However, in a common situation when most 
spins are frozen and Pfree ^ 1, one can neglect the dependence PfreeiT) and assume that 
1 — Pfree ~ 1- Wc argue, therefore, that the interference magnetoresistance mechanism is 
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not affected by small concentrations of impurities with free electron spin despite the fact 
that these impurities in the limit of low temperatures are included in the hopping process. 

Let us now discuss the saturation field for this negative magnetoresistance. The conven- 
tional theory which considers the destructive interference related to all tunneling paths leads 
to the saturation field Hgat ~ ^o/f'hP- It is the field that suppresses coherence in a charac- 
teristic pair of tunneling paths. Note that Hsat tends to zero with decreasing temperature. 
However, in our approximation for long hops one considers only significant paths that in 
some sense are similar to each other. Moreover, the interference phenomena are localized in 
the relatively small interference links. Accordingly, the saturation field is controlled by sup- 
pression of interference inside one link. It leads to Hgat ~ ^o/rscP- This field is independent 
on temperature and on hopping distance. 

There also exists the positive magnetoresistance due to spin alignment by the magnetic 
field, which is quadratic in H and has a strong temperature dependence. It arises from the 
term in (l3TI) related to the free spins and its expression is 

In ^ oc iVP,..,(ne^)/^a('^-)/^) (^)' oc r,P,..,T-i7^ (33) 

Again, this magnetoresistance mechanism has a strong temperature dependence. If the hop- 
ping is over the Coulomb gap states, the magnetoresistance is oc P/^ee^"^'^ where Pfree 
logarithmically decreases with temperature. The temperature dependence of this mecha- 
nism is stronger than the one for the other known positive magnetoresistance mechanisms 
including the wave-shrinkage magnetoresistance and the magnetoresistance due to the sup- 
pression of hops between the impurities with different occupation numbers 23|,|2J]. However, 
for the standard semiconductor with most electron spins frozen the value Pfree is small, so 
this magnetoresistance mechanism becomes important only at low temperatures. Also one 
should note that this mechanism saturates at magnetic field of the order of T/fibQ that is 
also small at low temperatures. 



A. numeric computation 

In our analytical theory we have used the approximation of large variance of scatterer 
energy. More exactly, only when the variance De of In Jq / | P | (where E is impurity energy 
measured from Fermi level) is large, the model of interference links becomes justified. Only 
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FIG. 7: (a) — interference contribution to magnetoresistance computed with eq. (112p for different 
numbers of scatterers N (x rh; (b) — dependance of linear expansion coefficient k on the scatterer 
number A^. The probabihty of an electron spin to be free was considered to be Pfree = 0.75. 

in this case one can discriminate between backbone, irrelevant and interference impurities 
and consider the number of interference impurities to be small. 
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FIG. 8: Interference contribution to magnetoresistance computed with eq. (|12p for different 
probabilities Pfree (a) and dependance of linear expansion coefficient k on 1 — Pfree (b). The 
calculation was made for resistors with = 10. 

In real semiconductors the discussed variance is of the order of unity. However, some of 
our results are expected to be valid for real semiconductors. For example, the log-normal 
distribution of tunneling path amplitudes remains correct for any De 0. Also we have 
shown numerically that some other results (the existence of significant pathes and backbone 
impurities) can be applied for realistic distribution of scatterer energies. Consequently we 
should also test our final results ( l32l) and ( l33l) with the numeric computations. 

In our computation we simulate the random distribution of scatterers in the area between 
the hopping impurities. We consider impurities only in the thin area around the hopping 
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line with the width p cx ^TscQ-. The energies of the scatterers are randomly distributed with 
constant density within the interval — Jq < Ej < Iq. Note that this distribution corresponds 
to rather small P^b. However we show that even for this distribution numerical results agree 
with our model. 

We consider the 2D distribution of the impurities. In our computation we keep the 
impurity concentration constant (na^ = 0.2), so that the number of scatterers is directly 
proportional to the hopping distance r/j. We assumed the hopping over the Coulomb gap 
states, so that ptg/T oc r^. 

To consider the system with free electron spins, we calculate all the tunneling path 
amplitudes Jj and estimate the magnetoresistance with a help of Eq. (IT^ . We performed 
the averaging over 1000 realizations (i.e. 1000 different critical hopping resistors). The 
results of these simulations are presented in figs. [7] and 13 We do not use any of our 
approximations (essential paths, backbone scatterers) in our computations. Instead we 
include all intermediate impurities and tunneling paths into the computation process. 

One can see that numerical results agree with analytic expressions (l32l) and ( l33l) . i.e. there 
is observable linear negative magnetoresistance (even for large Pfree)', the linear expansion 
coefficient k of magnetoresistance is proportional to cc N and 1 — Pfree- Also there is a 
positive quadratic magnetoresistance that increases with Pfree and increases with stronger 
than the linear magnetoresistance. 
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FIG. 9: The magnetoresistance for large hopping distances and Pfree = 0. (a) — magnetore- 
sistance for two different r/j. (b) — dependance of linear magnetoresistance on r/j and A^. One can 
see that dependance k oc rh stays justified when the number of scatterers grows faster than rh- 

However the question remains whether the computed linear dependence k{rh) is the real 
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semiconductor property or it is due to the model choice (in our computations we have chosen 
the scatterers in the strip of constant width, and the number of scatterers was oc rh)- 
We can not consider a larger number of scatterers and compute all tunneling amplitudes Jj 
because of the technical reason (the number of these amplitudes grow exponentially with 
A^). However we can bypass this problem when we consider the system without free spins. 
In this case we can use the summation algorithm that gives the correct result (in terms of 
expression (|T2|) ) with the computation time oc A^^. 

Accordingly we calculated the magnetoresistance for the case of Pfree = and large 
^h/O' ^ 35. We also took scatterers in a much wider area with the width p oc (so 
A^ oc rf/"^). This area contains all scatterers that can be included in tunneling path without 
making it exponentially small (independent on backbone impurities). Figure shows the 
corresponding results. One can see that the linear magnetoresistance is proportional to Vh 
and not to A^. Finally this calculation shows (fig. [9] (a)) no signs of decrease of saturation 
magnetic field with increasing hopping distance. This fact gives additional support to the 
model of interference links. 



VI. COMPARISON WITH EXPERIMENT 



In this section we discuss one experimental result that does not agree with the conven- 
tional picture of hopping negative magnetoresistance. These experiments demonstrate the 

3- 
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In the 



suppression of negative magnetoresistance at low temperature observed in 
corresponding experiments the magnetoresistance was measured at different temperatures 
in 2D GaAs-AlGaAs heterostructures, where both the wells and the barriers were doped 
by acceptor impurity Be and, thus, the acceptors within the well included double occupied 
A'^ centers. In contrast to the predictions of the conventional theory , it has been 

observed that the negative magnetoresistance is suppressed at low temperatures. 

It worth noting that the similar phenomenon has been observed earlier in bulk semicon- 
ductors [isl. It has been explained with the assumption of non-Coulomb character of 



impurity pot ential at large distances and effect of the Coulomb gap. However, as explicitly 
shown in [l9| , this explanation does is not valid for 2D systems (due to different asymptotic 
of the wavefunctions of localized electrons in 2D with respect to 3D one). So the spin align- 
ment magnetoresistance mechanism has been invoked to explain this phenomenon in 2D. 
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However no detailed theory of this magnetoresistance mechanism has been proposed. 




FIG. 10: Magnetoresistance of 2D semiconductor structures at different temperatures. The 
experiment results (dots) from [l^ are compared with parabolic law (I34p (lines). 

It is important to note that the systems observed in [l^ [3] include scatterers of two types 
(due to the double occupation) — the first type is the acceptor in the quantum well that has 
high activation energy and a short localization radius. The second type corresponds to the 
double occupied acceptors in the well 25|, |26| and the localized states of the hole in the well 
bound to the barrier acceptors. These states have relatively small activation energy and a 
large localization radius. It is likely that the spins of holes on single occupied acceptors in the 
well are free, while the spins of holes on the scatterers of the second type are frozen. Thus, 
in the considered system one has noticeable values of both Pfree and 1 — Pfree, and therefore 
both negative linear magnetoresistance and positive spin ordering magnetoresistance should 
be present in the system. 

Figure HO] shows these experimental data compared with the parabolic law 

which is simply the first two terms of the magnetoresistance expansion. One can see that 
the experimental data are in a good agreement with (1341) . It means that the higher terms 
of expansion can be neglected when discussing these experiments. 
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FIG. 11: The temperature dependence of k (a) and k2 (b) (dots) from compared to the 
theoretical results ([32]) and ([33|) (lines). 

The next step is to compare temperature dependencies of k and k2 with theoretical 
predictions. Assuming that the conduction is controlled by the Coulomb gap states the 
conventional theory predicts k oc T~~^/'^ and k2 oc T~^/'^ due to wavefunction shrinkage mag- 
netoresistance mechanism. These dependencies lead to the increase of the role of negative 
magnetoresistance at low temperatures. The extremal (minimal) magnetoresistance value 
increases with decreasing temperature as oc 

Our theory fl32]l33|) for long hops gives 
k oc T"^/^, also we predict the strong dependence ^2 oc T^^^^ due to spin alignment mag- 
netoresistance. The value of the magnetoresistance minimum depends on temperature as 

Figure [11] compares experimental data with our results fl32|33p demonstrating a good 
agreement with our predictions. The comparison with the conventional theory results in a 
significantly worse fit. So we argue that our theory is at least in semiquantitative agreement 
with experiment. 

VII. DISCUSSION 



As we have mentioned earlier, the conventional considerations of the interference mech- 



anism of magnetoresistance in the long hop limit 8l-ll2| neglect one important factor. It is 
the variance of logarithm of scatterer energy De- In the preceding studies this parameter 
has been set to be De = whether explicitly or implicitly. We started our theory from 
the opposite limit ^ 1. Actually, when De = oo, there is no interference magnetore- 
sistance as the single tunneling path dominates all the hopping process. Our theory can be 



27 



considered as a first correction related to a finite De- One can calculate further terms by 
considering sequences of two interference impurities separated by a backbone ones, then of 
three interference impurities and so on. 

We did not intend to go this far in the present study. However, we presume that the 
proportionality of magnetoresistance to should exist also for the higher terms. Actually, 
for any De ^ there is some probability to find a backbone impurity — the inclusion of 
this impurity significantly increase the tunneling amplitude. So the resistor in the long hop 
limit is still separated into the parts with this backbone impurities although these parts can 
be large for small (but finite) De- The magnetoresistance should be proportional to the 
number of these parts and thus should be proportional to r^. The dependence on Pfree can 
be different in the more rigorous theory as the parts separated by the backbone impurities 
can contain more than one free spin. 

One question is widely discussed in the physics of hopping with interference in the long 
hop limit. When the hops are short, the net tunneling amplitude is usually positive. It is 
not clear, however, whether the sign of net tunneling amplitude is completely random in 
the long hop limit or whether a positive sign still prevails. Our theory gives the following 
answer for this question. The sign is completely random when the resistor contains at least 
one backbone impurity. Naturally the net tunneling amplitude is controlled by significant 
tunneling paths that include this backbone impurity. The sign of the energy of this impurity 
is random, so the sign of the net tunneling amplitude is also completely random. Clearly, 
at any temperature (and hopping distance rh) there is some chance to find a resistor with- 
out backbone impurities but the probability of such an event drops exponentially with rh- 
Therefore the average sign of net tunneling amplitude also drops with Vh exponentially. 

Summing up, we generalized the theory of interference effects in the hopping magnetore- 
sistance to include the contribution of scatterers with free electron spins. We considered 
both the case of short hops when the mean number of scatterers in a resistor is less then 
unity (many resistors in this case have no scatterers and thus no interference effects) and 
the case of long hops when the mean number of scatterers in a resistor is large. For the case 
of long hops we developed a new approach to the problem of interference magnetoresistance 
that is based on assumption of large variance of the logarithm of scatter energy De (the 
conventional approaches consider De — 0). We showed that our approach is in a good 
agreement with numerical computations with realistic value oi De- Our theory allows one 
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to calculate explicitly the temperature dependence of magnetoresistance and its dependence 
on the relative ratio of free spins Pfree for the short hops and the long hops cases. Our 
results are in semiquantitative agreement with experimental data on magnetoresistance in 
GaAs — AlGaAs 2D structures. 

We are grateful to A.S. loselevich and N.S. Averkiev for many fruitful discussions. Our 
work was supported by RFBR foundation (project n. 10-02-00544). A.V.S. also acknowl- 
edges support from "dynasty" foundation and from the program of RAS presidium "the 
support for young scientists". 



Appendix A: The expression for with the assumption of large number of back- 
bone scatterers and with an account of free spins. 



To consider the problem of magnetoresistance in a system with free spins, we should start 
from the expression (fTOl) for the hopping rate. We also adopt the model of independent links, 
i.e. we assume that is controlled by the significant tunneling paths that differ from each 
other only in local interference links. 

First, we note that inclusion or exclusion of a scatterer with frozen spin does not changes 
the tunneling result independently of spin configuration. Each significant tunneling path j 
has a pair that differs from j only by the interference link i (for any link i). One can see 
that if link i does not contain a free spin, then by combining such pairs one can take out 
the factor 



Iil,i3 + 

corresponding to this interference link i. 



Iil'i2li2. 



i3 



i2 



(Al) 



Also one can take out all factors corresponding to backbone scatterers. However, these 
factors do not contribute to interference magnetoresistance, so we omit them 
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Here the first product is taken over interference links with frozen spin or interference links 
with unoccupied impurities. The product over s is taken over interference scatterers with 
free spins. L^J^ is the amplitude of the interference link s in the path j. It can be equal 
to either /si,s3 if the path j does not include scatterer s or to Isi,s2ls2,s3/Es2 if the path j 
includes s. 
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Further, we use the approximation A*";,;, 3> A^mt once again. We consider that for any 
interference hnk s with free spin the previous impurity with free spin s— is a backbone one. 
Accordingly the next impurity with free spin s+ (for any interference hnk with free spin s) 
is also a backbone impurity. 

With this assumption one can see that two paths that differ only by a single impurity 
with free spin s lead to the same tunneling result for all configurations when spin projections 
on s and s— are the same. In other spin configurations this tunneling paths lead to different 
results, however these results differ only by the replacement of final projections of spin on 
the impurities s and s+. Note that in this model the inclusion or exclusion of an interference 
impurity s does not change the interference pattern on the other interference impurities with 
free spin that follow s. This interference is controlled by initial spin projections on s+ and 
on other backbone impurities with free spins. Thus one can obtain the following expression 

for ri2. 
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Here P^^ is the probability for spins on impurities s and s— to have the same projection. 

= Pp^ with Pp defined in (JT]). 



Naturally, this probability does not depend on s: P, 



tt 



Similarly, P. 



n 



1 — Pf^ for any interference link with free spin s. 
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